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ABSTRACT 

A statistical method is developed to infer the typical density profiles of poor 
galaxy systems without resort to binning of data or assuming a given center to each 
system. The method is applied to the accordant redshift quartets in Hickson's compact 
groups (HCGs). The distribution of separations in these groups is consistent with a 
unique generalized modified Hubble surface density profile, with best-fit asymptotic 
slope (3 = 1.4 and core radius Rc = 18/i~^kpc, although a King profile {(3 = 1) is 
also consistent with the data (with Rc = 6/i~^kpc). These distributions are more 
concentrated than what has been previously determined for these groups. HCGs are 
unlikely to act as strong gravitational lenses, but analogous systems 5 to 10 times more 
distant should produce a non-negligible fraction of giant arcs. 
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1. INTRODUCTION 

The 100 compact groups (CGs) cataloged by Hickson (1982) have been among the most 
systematically observed galaxy systems, and yet, there is little consensus on the nature of these 
systems, between those who believe that most of Hickson's compact groups (HCGs) are roughly as 
dense in 3D as they appear to be in projection (e.g. Hickson &; Rood 1988), and those who argue 
that CGs are mostly chance alignments of galaxies along the line of sight, within larger 'loose' 
groups (Mamon 1986) or cosmological filaments (Hernquist, Katz h Weinberg 1995). 

The nature of CGs may be clarified through the characterization of their galaxy surface 
number density profiles, S(i?). Indeed, if CGs are caused by chance alignments within much 
larger systems, then the galaxy positions should be relatively random, and the corresponding 
density profiles should be roughly homogeneous. The galaxy number density profile in CGs is 
also interesting by itself, as compared with the density profiles of the gaseous and dark matter 
components, in connection with the theory of the formation and evolution of galaxies within 
groups. 

The small number of galaxies (typically 4) per HCG makes it very difficult to estimate 
the density profile. Hickson et al. (1984) have shown that all 100 HCGs, rescaled to the same 
size, show a concentrated density profile. More recently, Mendes de Oliveira & Giraud (1994, 
MG94) have argued that HCGs, once rescaled to a common size, are consistent with a King 
(1962) surface density profile, with typical core radius of 15^~^kpc. Unfortunately, such studies 
suffer from the uncertainty introduced by the definition of a center for a system of typically four 
galaxies, and from the assumption that this center is at the galaxy barycenter. Nevertheless, 
there are statistical methods to determine the density profiles of discrete systems that do not 
make use of the definition of the system center, relying instead upon the distribution of projected 
particle (galaxy) separations (see the mathematical foundations in Salvador-Sole, Sanroma &: 
Gonzalez-Casado 1993, SSG; Salvador-Sole, Gonzalez-Casado & Solanes 1993, SGS). 

Alternatively, one can obtain good estimates of the surface density profiles of the X-ray 
emitting diffuse intergalactic gas, which, so far, has been firmly detected with the ROSAT PSPC 
in 7 HCGs out of 17 pointings (see Mulchaey et al. 1996 [MDMB], and references therein), 
suggesting that those HCGs are physical bound entities. These ROSAT observations reveal 
projected gas density profiles consistent with Sgas(-R) = Sgas/ll + {R/ RcY]^^^^'^''~^^^'^ , and with 
core radius 4 < i?c < 30 /i^^ kpc and 0.38 < /3gas < 0.92 (see MDMB, and references therein). In 
comparison, Abell clusters have X-ray core radii typically 15 times larger (Jones & Forman 1984; 
Mushotzky 1994). 

The efficiency of a structure as a gravitational lens is determined by its projected mass 
density profile, Emass(-R)) and the relative lens-observer and source-observer distances (Bourassa & 
Kantowski 1975). Extrapolating the small intergalactic gas core radii to the distribution of total 
binding mass suggests that CGs may be concentrated enough to be candidates for gravitational 
lenses (MG94). If the binding mass in groups is more concentrated than the intergalactic gas, as 
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is the case in clusters (Henriksen 1985; Hughes 1989; Durret et al. 1994), the lensing efficiencies of 
groups are enhanced. 

Hence, various signatures of gravitational lensing (giant arcs, arclets, shear, apparent 
magnitude enhancement) in CGs of existing catalogs (Hickson et al. 1992; Prandoni, lovino & 
MacGillivray 1994) or of future deeper surveys, if observed, would be of major interest because 
they would allow a quantification of the dark matter distribution in CGs through a method 
independent of X-ray analyses. This quantification has important consequences for the theories 
of formation and evolution of galaxies and their groupings in the framework of the large scale 
structure of the Universe. 

In this Letter, we implement a statistical method to determine average galaxy density profiles 
in the case of poor statistics, with no assumed group center. The mean galaxy surface number 
density profile of compact groups is computed from these statistical methods, and the implications 
for the lensing properties of compact groups are studied. 



2. THE GALAXY SURFACE DENSITY PROFILE 

For a catalog of groups (whose centers are uncorrelated), the galaxy two-point correlation 
function, ^(s), is (Peebles 1980): 

(iV/(S,.S,)(.))nG 

= :^ ' (1) 

where Sj(s) is the projected galaxy density profile of the j— group normalized to unity, s is the 
projected intergalactic separation, Ni is the number of galaxies in the i— group, no and n are the 
mean number density of groups and galaxies respectively, ( ) is the average taken over the groups 
of the catalog and (Sj * Sj)(s) is the convolution of the projected density profile of the 2— group, 
which can be calculated by means of the deconvolution method (SSG; SGS): 

(S*S)(s)=.Fio^|^o^fi[2P(>s)]} , (2) 

where J-i and A are the ID Fourier and Abel transforms, respectively, and P{> s) is the fraction 
of galaxy pairs with intergalactic separation > s. 

The average galaxy surface density profile for groups in a given catalog can be determined 
through a fitting procedure by means of equations (||) and (§) above. Assume that the groups in 
the catalog are discrete realizations of an analytical profile, S(i?). We use here the generalized 
modified Hubble profiles (hereafter, GMHP, see Jones & Forman 1984), whose shape is characterized 
by the core radius, Rc, and the slope, /?: 



(3) 
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where Sq is the central projected density. Our free parameters Rc and P are obtained by 
calculating ^^ata^^-j £qj, ^.j^g data catalog and ^'^ {s; Rc, for the analytical profile, using equations 
(|l]) and (^), and minimizing 



where A'sep is the number of bins of intergalactic separations, s^- The dispersions = ((Cfe — (Cfc))^)) 
required in evaluating equation are evaluated from Monte-Carlo simulations of artificial group 
catalogs, generated with GMHPs, with the same number of galaxies per group and of groups per 
catalog as the real group catalog. Since the different are not independent random variables, 
the number of degrees of freedom of X'^ is not known, and we have to resort to Monte-Carlo 
simulations to evaluate empirically its distribution (instead of the usual statistics), for 
estimating confidence levels. Note that only information on the projected intergalactic separations 
is required in the deconvolution method. This circumvents the delicate problem of determining 
the center of the groups. Moreover, noise suppression by filtering is also implicit in this method. 

In catalogs of groups, one has a low number of galaxies per group and of groups per catalog. 
In this case of poor statistics, the reliability of the method has been tested with simulations of 
artificial catalogs whose groups are realizations of GMHPs with fixed /?. Our results show that the 
method is able to recover the input value of Rc to a typical 90% confidence precision of a factor 
1.6, and to distinguish between different values of (3 in the case of poor statistics (Montoya et al. 
1995). 

The method explained above has been applied to the HCG sample of 42 accordant redshift 
quartets (Hickson et al. 1992). We have tested the null hypothesis Tio that the groups of this 
sample are realizations of a unique GMHP with 1 < /3 < 3. For each value of (3, we have calculated 
X"^ (eq. 0) for the sample, and the cumulative probability F(X^) from Monte-Carlo simulations. 
The resulting confidence levels in the {Rc, (3) plane are plotted in Figure 1. 

The optimum fit is obtained for fi = 1.4 and Rc = 18 kpc (hereafter, model 2), 
corresponding to F(X^) = 0.16. For /3 = 1,2, and 3, the best fit is for Rc = 5.7,32.7, and 
50.1 kpc (models 1, 3, and 4 respectively). Our fitted surface density profiles have been 
obtained without any rescaling of the sizes of the groups. Such rescalings are not necessary 
because we get good fits without them. To show this, we have compared the distributions of the 
minimum, median and maximum intergalactic separations for each group in the catalog, with 
those for artificial catalogs of 42 quartets, each of them being a realization of a unique GMHP 
with P = 1.4 and Rc = 18 kpc. Table 1 shows the median Kolmogorov-Smirnov probabilities 
of not being able to reject the hypothesis that the distribution of separations of the catalog are 
consistent with that of the Monte-Carlo catalogs with (3 = 1.4 and Rc = 18 kpc. We conclude 
that any rescaling (using the minimum, median and maximum intergalactic separations as typical 
lengths) leads, indeed, to much poorer fits, whatever values of both /5 and Rc- 



MG94 have also computed the galaxy surface number density profile for a sample of HCGs. 
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Their fitting procedure requires a group center (for which they adopt the barycenter), and, 
moreover, they scale all groups in their sample to a common size by normalizing each group to 
the largest galaxy distance to the group barycenter, i?max- Both of these operations introduce 
important uncertainties in their results. On one hand, the center of a system of few (typically 4) 
galaxies is not well defined. Beers and Tonry (1986) have shown that any possible inner galaxy 
concentration is washed out when the barycenter is used as the center of galaxy clusters. So, given 
P = 1, MG94's core radius should be much larger than ours, and indeed it is 2.5 times larger (Fig. 
1). On the other hand, it is not clear that groups should have a common size and that i?max 
should be a useful measure of that size. Instead, our results show that our sample of HCGs can 
be considered as a realization of a unique profile, and that the dispersion in .Rmax arises from 
statistical fluctuations. 

The galaxy density profiles found here are steeper than the gas density profiles obtained from 
X-ray surface brightness data (see § 1). Qualitatively, this can be expected in the framework of the 
/3-model (Cavaliere 1974), which predicts Pspec = <''u/(^^/a*'^p) = Pgas/Pgal- Indeed, the galaxy 
system is cooler than the X-ray emitting gas, as HCGs have /3spec — 0-7 (Mamon & Henriksen 
1996). The typical values of /3gas = 0.6 and our best fit of /3gai = 1.4 are consistent with /3spcc in 
the context of the /3-model, given the uncertainties in these three parameters. Hydro dynamical 
simulations of X-ray cluster formation including galaxy formation (e.g. Frenk et al. 1996) also 
predict a velocity bias and spatial gas/galaxy segregation. 



3. LENSING BY COMPACT GROUPS 

Due to their peaked density distributions (§ 2), HCGs could constitute efficient gravitational 
lenses. An HCG at angular distance Di with central projected mass density Smass(O) may produce 
an elongated arc image of background sources at a distance Dg if (Bourassa & Kantowski 1975) 
Smass(O) > (c^/47rG) {Dg/ [DiDis]) (strong lensing condition, hereafter SLC). For Q = 1 and a 
smooth universe, the angular distance is Di = 2{c/Ho){Gi — 1)/Gf , where Gj = (1 + Zi)^/"^, i = l,s, 
and one can write (Blandford & Kochanek 1987) A, = 2{c/Hq){Gs - Gi)/{GlGi). The SLC 
imposes, for a given lens, a minimum source redshift: 



mm 



AGi{Gi-l)-Gt 
A{Gi-l)-Gf 



,4n2 



-1, (5) 



where A = 87rGSmass(0)/(c-ffo)- A. physical solution z™"^ > zi is obtained, after some algebra, 
directly from the SLC when Smass 

(0) > (1 + ziY/'^ [c^/ {AttGDi)]. Assuming that the galaxy 
system has isothermal and isotropic kinematics, one can solve for the total mass density in Jeans 
equation, and integrating over the on-axis line-of-sight, one finds Sinass(O) = 3f3(j'^/{2GRc), where 
ay is the group ID velocity dispersion. Dropping the factor (1 -|- zi)^/"^ ~ 1 -|- zi/2, the condition 
for zf^^^ > zi can now be written: 

2 R. 
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In Figure 2, we plot {thick histograms) the cumulative fraction of HCGsQ from our sample 
for which the product Di{ay/c)'^ exceeds different thresholds (eq. [^), marked by arrows. For 
model 1, three HCGs (8, 60, and 82) in our sample verify equation (^), implying (eq. [|5|) 
Zs > 0.19,0.20 and 0.08, respectively. For our other models, no HCG in our sample fulfills the 
strong lensing condition. The three HCGs that are potential strong gravitational lenses have 
much higher line-of-sight velocity dispersion than the average (517, 472, and 714kms~^, compared 
to a median of 237kms~^ for our HCG sample), and this could easily be due to sampling and 
projection effects, as is suggested by the lack of HCGs detected in X-rays with gas hotter than 
IkeV (MDMB). Put another way, since most HCGs detected in X rays have temperatures close 
to 0.9 keV (see MDMB), with Mamon & Henriksen's (1996) estimate of /^speo we infer true ID 
velocity dispersions of = 315kms~^. The thin histograms in Figure 2 illustrate the situation if 
all HCGs have ay = 315 km s^^, and we infer that no HCGs induce strong gravitational lensing. 

HCGs are relatively nearby (^median = 0.03) and future catalogs should provide more distant 
CGs. In light of this, we have considered a sample of CGs with the same properties as before, but 
with Dia'^ multiplied by a factor of 5 and 10, and we have plotted in Figure 2 the corresponding 
cumulative fractions of CGs whose Di{ay/cY values exceed the same thresholds as before. For 
model 1, the strong lensing criterion is fulfilled for 38% and 62% of equivalent CG samples 5 
and 10 times more distant, respectively, using the measured velocity dispersions for a^, and 69% 
and 95%, respectively, using the canonical value = 315kms~^. Thus, we obtain more efficient 
lensing in comparison with MG94, who, with the same profile, predict that 37% of a sample of 
HCGs 10 times more distant will produce strong lensing. For our best fit model 2, we expect 
strong lensing for 18% and 38% of our sample, respectively (see Fig. 2), using the measured 
velocity dispersions, and 21% and 64% with = 315kms~^. These results confirm MG94's 
suggestion that CGs in deeper surveys (and/or with larger velocity dispersion than for HCGs) 
could turn out to be efficient strong lenses. 

Lens systems cause an enhancement of the apparent magnitude of a source placed behind 
them. Assuming a perfect alignment and Dg = 2 Dig, the fraction of HCGs that amplify 
background sources by > 0.5 mag is 21% for model 1 and 8% for models 2 and 3, using measured 
UyS, and 26% for model 1 and negligible for the others, using = 315kms~^. In comparison, 
MG94 find 16%. 

4. DISCUSSION 

The uniqueness of the absolute surface density profiles of HCGs may seem puzzling when 
considered within the hierarchy of galaxy systems. It arises from the characteristic mass density 
and velocity dispersion required for virialized groups from a combination of sharply decreasing 
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cosmological distributions of these properties and minimum cutoffs for virialization (Mamon 1994). 
Consequently, virialized groups have a characteristic radius, and by extension, a characteristic 
density profile. 

Because our method for determining the density profiles of groups does not make use 
of a choice for a group center, only relying on the projected separations, our fitted surface 
density profiles are more concentrated than that obtained by MG94. Our profiles are even more 
concentrated than for most HCG diffuse intergalactic gas distributions, as determined from X-rays. 
Part of this inferred concentration may be due to small-scale correlations, i.e., binarity, in the 
galaxy distribution (see Walke & Mamon 1989). Note that, whereas Hickson (1982) only selected 
the bright galaxies, the spectroscopically confirmed faint members of HCGs have somewhat larger 
projected separations (Zepf, private communication). 

Although the lensing efficiency of HCGs is highly model dependent, we confirm MG94's result 
that HCGs should not constitute strong lenses, as confirmed by the absence of giant arcs in CCD 
images of HCGs (see Hickson 1993). CGs in 5 to 10 times deeper surveys should produce a few 
giant arcs. Our results do not exclude that a fraction of HCGs are due to projection effects from 
loose groups or long filaments. In this case, real groups in Hickson's catalog should be even more 
concentrated, and our conclusions on lensing would be reinforced. 
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Table 1. Kolmogorov-Smirnov test probabilities 
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Fig. 1. — Best-fit surface density profile (eq. [3]) parameters of HCGs, for eonfidence levels of 
100%, 95%, 90% and 85% {solid curves) and curve of optimum fits {dashed curve). Numbers refer 
to our adopted models (with '2' our best solution). 
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Fig. 2. — Cumulative fraction of HCGs with values of Di(ay/c)'^ greater than given thresholds (eq. 
[6]), for our sample of HCGs, and for the same samples placed at 5 and 10 times their original 
distances. Arrows show the thresholds for strong Icnsing satisfied for models 1 {dashed arrow), 2 
{solid arrow), and 3 {dotted arrow). The threshold for model 4 is similar to that for model 3. Thick 
histograms use measured velocity dispersions, while thin histograms use o"„ = 315kms~^. 



